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Abstract. Cosmic rays accelerated by a shock form a streaming distribution 
of outgoing particles in the foreshock region. If the ambient fields are negligible 
compared to the shock and cosmic ray energetics, a stronger magnetic field can 
be generated in the shock upstream via the streaming (Weibel-type) instability. 
Here we develop a self-similar model of the foreshock region and calculate its 
structure, e.g., the magnetic field strength, its coherence scale, etc., as a function 
of the distance from the shock. Our model indicates that the entire foreshock 
region of thickness ~ ij/(2r^j^), being comparable to the shock radius in the late 
afterglow phase when Fsij ~ 1, can be populated with large-scale and rather strong 
magnetic fields (of sub-gauss strengths with the coherence length of order 10^^ cm) 
compared to the typical interstellar medium magnetic fields. The presence of such 
fields in the foreshock region is important for high efficiency of Fermi acceleration 
at the shock. Radiation from accelerated electrons in the foreshock fields can 
constitute a separate emission region radiating in the UV/optical through radio 
band, depending on time and shock parameters. We also speculate that these 
fields being eventually transported into the shock downstream can greatly increase 
radiative efficiency of a gamma-ray burst afterglow shock. 
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1. Introduction 



Do gamma-ray bursts (GRBs) accelerate cosmic rays (CRs)? There are arguments 
in favor of such an idea (Dermer & Atoyan 20041. It is supposed that CRs are 



accelerated via the Fermi mechanism in which a particle crosses the shock many times 
and gradually gains its energy. For an ultra-relativistic shock, a steady-state universal 
power-law energy distribution of particles shall form ( Kirk et ar]|2000 Achterberg et 
al.|2001 1. The shock shall accelerate all particle species, but the electrons being much 
lighter than protons and ions will also loose their energy via synchrotron cooling. This 
radiation is thought to be observed as the delayed afterglow emission of GRB sources. 
A problem immediately arises here from a simple estimate: if the X-ray afterglow 
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observed on a day timescale after the prompt burst is indeed the synchrotron radiation 
from the shock-accelerated electrons, then the pre-shock medium has to be highly 
magnetized with the fields of milligauss strengths ( Li fc Waxman|[2006 1 . Thus, either 
the magnetic field is somehow generated in the shock upstream, or the conventional 
paradigm of the GRB afterglow needs revision. 

In this paper, we present a self-similar model of the large region in front of a 
relativistic shock - the foreshock. This region is populated with the shock-accelerated 
particles, which stream away from the shock into the collisionless ambient medium and 
generate magnetic field via a streaming (Weibel-type) instability. The model predicts 
the generation of strong, sub-gauss, magnetic fields in the entire foreshock whose 
thickness is ~ R/{2V1y^ and is comparable to the shock radius, ^ 10^'' — 10^® cm, in the 
afterglow phase a day or more after the explosion, when the shock is weakly relativistic 
or non- relativistic. The fields are sustained against dissipation by the anisotropy of 
newly accelerated particles. Moreover, these fields are relatively large-scale, with the 
coherence length being as large as a fraction of the foreshock size, ~ 10^^ cm, which 
makes them effectively decoupled from dissipation. We speculate, that these mesoscale 
magnetic fields being ultimately advected into the shock downstream can significantly 
increase the radiative efficiency of GRB afterglows and, perhaps, explain the origin 
of the magnetic field in an external shock of a GRB. We remark here, however, that 
our study is analytical and cannot account for a number of nonlinear feedback effects 
of the generated fields and pre-conditioned external medium onto the shock structure 
and particle acceleration. Kinetic and hybrid computer modeling is essential for better 
and more accurate understanding of the foreshock structure. 



2. The model 



Overall, our model is as follows. A shock is a source of CRs which move away 
from it, thus forming a stream of particles through the ambient medium, say, the 
interstellar medium (ISM). If the ISM magnetic fields are negligible, i.e., their energy 
density is small compared to that of CRs, the streaming instability (either the 
pure magnetostatic Weibel or the mixed-mode electromagnetic oblique Weibel-type 
instability, depending on conditions) is excited and stronger magnetic fields are quickly 
generated. These fields further isotropize (thermalize) the CR stream. Since less 
energetic particles, having a greater number density and carrying more energy overall, 
are thermalized closer to the shock, the generated B-field will be stronger closer to the 
shock and fall off away from it, whereas its correlation length will increase with the 
increasing distance from the shock. More energetic particles keep streaming because of 
their larger Larmor radii and produce the magnetic field further away from the shock. 
This process stops at distances where either the CR flux starts to decrease (because of 
the finite distance the CR particles can get away from a relativistic shock or because 
of the shock curvature causing CR density to decrease as oc if the shock is sub- 
or non-relativistic) or where the generated magnetic fields become comparable to the 
ISM field and the instability ceases. Thus, a large upstream region — the foreshock 
— is populated with magnetic fields. We now derive its self-similar structure. We 
work in the shock co-moving frame unless stated otherwise. 

Let's consider a relativistic shock moving along a;-direction with the bulk Lorentz 
factor Fsh; the shock is plane-parallel and lies in the yz-plane, and a; = denotes 
the shock position. The shock continuously accelerates cosmic rays, which then 
propagate away from it into the upstream region. We conventionally assume that 
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the CR distribution over the particle Lorentz factor is described by a power-law: 

"•CR = ?^o(7/7o)"'' (1) 

for 7 > 7o and zero otherwise. Here the index s = p — 1 is approximately equal 
to 1.2 for ultrarelativistic shocks and no is the normalizationj|] We assume that the 
above energy distribution is the same everywhere in the upstream, that is, we neglect 
the nonlinear feedback of magnetic fields onto the particle distribution. The CR 
momentum distribution exhibits strong anisotropy: the parallel {x) components of 
CR momenta are much greater than their thermal spread in the perpendicular {yz) 
plane. Indeed, for a particle to move away from the shock, it should have the a;- 
component of the velocity exceeding the shock velocity. Since both the shock and 
the particle move nearly at the speed of light, this puts a constraint on their relative 
angle of propagation to be less then l/Fsh in the lab (observer) frame. Hence, the 
transverse spread of the CR particle's momenta is p±^ < p\\ /T^h ^ V\\ ■ This is also 
seen in numerical simulations ( Spitkovsky||2008 ) . 



The CR particles propagate through the self-generated foreshock fields and scatter 
off them. Lower energy particles are deflected in the fields more strongly and, 
therefore, izotropize faster than the higher energy ones, as having larger Larmor 
radii. At a position a; > the CR distribution can roughly be divided into isotropic 
(themalized) component with 7 < 7^(2^) and streaming component with 7 > ^r{x), 
where ^r{x) is the minimum Lorentz factor of the streaming particles at a location 
X] it is also the maximum Lorentz factor of the randomized component at this 
location. The streaming component is Weibel-unstable with a very short e-folding time 

1 /2 

T ^plei, where Wp^rei = (47re^fi(7)/mp7) is the relativistic plasma frequency, h{"f) 
is the density of streaming particles of the Lorentz factor 7 (tilde denotes streaming 
particles). Note that the Weibel instability growth rate depends on n of the lower 
density component - cosmic rays, in our case - measured in the center of mass frame 
of the streaming plasmas. For the lower-energy part of the CR distribution, the 
center of mass frame is approximately the shock co-moving frame, hence we evaluate 
the instability on the shock frame. This approximation is less accurate for the high- 
energy CR tail; however, the growth rate and the scale length are weak functions of 
the the shock Lorentz factor (cx L^^^^), so the result will be accurate within an order 
of magnitude for all reasonable values of V^h for GRB afterglows. Here we use the 
proton plasma frequency because the CR electron Lorentz factors are about rrip/me 
times larger, so they behave almost like protons ( Spitkovsky||2008 ) . The instability is 



very fast: it rapidly saturates (the fields cease to grow) in a few tens of e-folding times 
r, that is in few tens of inertial lengths (also referred to as the ion skin length) c/wp_rci 
in front of the shock. Thereafter the particles keep streaming in current filaments and 
the field around them amounts to ^ 0.01 — 0.001 or so of the kinetic energy of this 
group of particles: 

^'(7)/87r^^smpcV(7), (2) 
where is the efficiency factor obtained from particle-in-cell (PIC) simulations; 
has the same meaning as the conventional ts parameter reserved here for the ratio 
of the total magnetic energy to the total kinetic energy of the shock and which, as 
is seen in PIC simulations, is larger than near the shock because of the nonlinear 

§ Conventionally the distribution is given as dn/d'y oc with p being ~ 2.2 — 2.3 for relativistic 
shocks; hence the density of particles of energy ~ 7 is n(7) oc ■y~PS'y oc ■y~P'^^ . 
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evolution and filament mergers. The correlation length of the field is of the order of 
the ion inertial length 

A(7) ~ c/ujp^rci = (mpC^7/47re^n(7))^^^ . (3) 

These random fields deflect CR particles and ultimately lead to their 
isotropization. The deflection angle of the particle on a field coherence length scale in 
the self-generated field is 

0^Sp^/p^iX/c)u;B, (4) 

where lob — eB/jmpC, p is the particle momentum and 5p± is it's transverse change. 
Using Eqs. ([s]), we obtain: 

e ~ eB(7)A(7)/(7"^pc2) ^ ^/2^. (5) 

Note that the deflection angle is independent of the particle's energy, as long as the 
field is produced by the particles of the same energy 7. The particles diffuse in the fleld 
and their rms deflection angle after transiting through a distance a; is ^ Q^fxJX. 
The group of particles thermalizes when '--^ 1, i.e., at the distance from the shock: 

Xr ^ X/e^ ^ A(7)/(2Cb). (6) 

At this point, x — Xr, one has 7 = 7r by definition; no field of the strength -B(7r) 
and the scale A(7r) can be produced at a; > Xr- Similarly, one can estimate the 
randomization of the higher energy particles with 7 3> ^r'- ^(7) ~ V^Sjsilr/^) <C 
^(7r)j which means that these particles keep streaming through much larger distances 
X ^ Xr and will produce the magnetic field further away from the shock. This field 
will be weaker and larger scale because of the lower density of the streaming particles 
n(7) n{'-fr), according to Eqs. ([T])-(|3|. 

Finally, the number density of streaming CR particles at 7,. is ^(7^) = 
"-o(7r/7o)"''- Therefore, 

A(7.) ~ (mpC^o/Wno)'/' (7./7o)^'+^^/' ^ Ao(7./7o)<'+^^/', (7) 

where Aq is the inertial length of the lowest energy CR "plasma". Inverting this 
expression yields: 

Ir - 7o[A(7.)/Ao]'/(^+^) ~ 7o(2^Ba:./Ao)2/(i+^). (8) 

Hereafter, the subscript "r" can be omitted without loss of clarity. 

In a steady state, this field is continuously advected toward the shock (in the 
shock co-moving frame since the center of mass frame of the foreshock plasma differs 
from the shock frame) and may affect the onset and the saturation level of the Weibel 
instability. In addition, the current filaments producing the fields merge with time, 
so that B and A change while being advected. These nonlinear feed-back effects are 
difficult to properly account for in a theoretical model; hence they are omitted in the 
current study. PIC simulations can help us to quantify the effects as well as to confirm 
or disprove our assumption that the shock and the foreshock do form a self-sustained, 
steady state structure. 
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3. The self-similar foreshock 

The self-similar structure of the foreshock immediately follows from Eqs. ([l]), ([2|, 
^ and Q. The magnetic field correlation length is proportional to the upstream 
distance from the shock, 

A(x)~x(2eB), (9) 
and its strength decreases with the distance as 

B{x) ^ Bo{x/xoy^ , (10) 

where Bq = {STr^BrripC^nojo)^^^ and a;o = Xo/{2^b) = {nipc'^^o/^-^e'^noy^^ /{2^b)- 
In this estimate we neglected the advected fields -6(7) as sub-dominant compared 
to B{'-fr) for -f > jr- The parameter expresses the field energy normalized to the 
shock kinetic energy. The energy of cosmic rays is Uqr = J / ^o){mp(?^) ^(7/70) ~ 
rripC^ jqUq and constitutes a fraction ^cR of the total shock energy, Ush- The efficiency 
of cosmic ray acceleration, ^cR, can be as high as several tens percent, perhaps, up 



to ^CR ~ 0-5, as follows from the nonlinear shock modeling (Vladimirov et al. 2006 
Ellison et ar]|2007| . The scaling of is: 



es -CcrCb (x/xo) . (11) 

These scalings hold while the shock can be treated as planar and while the ISM 
magnetic fields are negligible compared to the Weibel-generated fields. If the shock 
is relativistic, CR particles can occupy a narrow region in front of it. Assuming CR 
to propagate nearly at the speed of light, their front is ahead of she shock at the 
distance Ar — ct^-d = c{R/c — R/vsh) — R — R/[i^ — ^f^^lh)] measured in the lab 
(observer) frame, that is at the distance ^ Ar/Fgh ~ ^/(2rsh) in the shock frame. 
Also, when the radial distance in the lab frame Ar = x/Tgh becomes comparable 
to the shock radius Ar ^ R the curvature of the shock can no longer be neglected: 
the density of CR particles, which was assumed to be constant in our model, starts 
to fall as cx r~^. This leads to a steeper decline of B with distance. Obviously, 
the first constraint is more stringent for a relativistic shock, whereas both are very 
similar (within a factor of two) for a non-relativistic shock. Hence we use the first 
constraint hereafter. Meanwhile, at some distance X, the Weibel-generated fields can 
become comparable to the ambient magnetic field, B{X) ~ -Bamb and the Weibel 
instability ceases; here we used that the ambient field in the shock frame is i?amb ~ 



-BisMj-LTsh ~ SisMrsh- PIC simulations (Spitkovsky 20051 indicate that for low 



magnetizations a < 0.01, i.e., _B(X)/i3amb > 0.1, the shock behaves as unmagnetized 
and the Weibel instability dominates. Although there is no a sharp threshold, one 
sees the Weibel instability to be suppressed for lower values of _B(X)/i3amb- To 
the order of magnitude, we set BisuTsh - B{X) - Bo{X/xo)''^'-^^/'-'+^\ therefore 
X - xo [Bo/{BisMTsh)]^'^^^^^'~^^ ■ To conclude, the the scalings, Eqs. ([9])-([ll}, hold 
at a; < Xmax, where 

x„,ax - Min [i?/(2r3h), X] - Min [i?/(2r,h), xo {Bo/ BisMTsh)"^] -(12) 

The region filed with the magnetic field in front of the shock is large, so does 
the region where radiation emitted by the CR electrons. The power emitted by a 
relativistic electron in a magnetic field is Pb = (4/3)crTC7e (^^/Stt), where ax is the 
Thompson cross-section and 7e is the Lorentz factor of the emitting electron. This 
expression is accurate for both synchrotron and jitter radiation ( Medvedev|2000 1 . For 
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the distribution of electrons (jlj homogeneously populating the foreshock, the power is 
dominated by the lowest energy particles with 7e ie{'mp/'>TT'e)jp ~ ee(™p/we)7o ''^ 
ee{mp/mi,)Tsh (here Ce is the efficiency of the electron heating) and the density 
Ue '-^ rip ^ riQ ^ nisMTsh- Then Ptot = / PBix)nodV , where the volume element 
is dV — inR'^dx, so 

/■^"•■■'- _^p-l 3^ 

Ptot OC / {x/xo) P+'-dx ^ Xo{Xynax/xo)^ + ' ■ (13) 

Jo 

Note that the co-moving radiated power increases with the foreshock thickness, 
Ptot oc a^max'''^'"^^'''', thus cmissiou is not localized in a thin layer near the shock 
and is, in fact, dominated by large distances: 

Ptot ~ LcR^Ble i2/3)aTnoXQ {Xn^ax/Xo) , (14) 

where Lqr = 47rP^(rnpC^no7o)c is the kinetic luminosity of cosmic rays, which is 
a fraction ^cR < 1 of the total kinetic luminosity of a GRB, — ^cr^grb, R 

is the shock radius and the co-moving CR density is of order the density of the 
incoming ISM plasma, no ~ "-igMrsh- The foreshock electrons are radiating in 



the synchrotron regime: the jitter parameter S (Medvedev 2000 Medvedev et al. 



2007|, which is the average deflection angle of an electron in the foreshock fields, 

6e ~ {eB{x)/"femec){X{x)/c) over the radiation beaming angle, ~ l/7e, is always 
much larger than unity: 

S{x) ^ eS(x)A(x)/(m,c2) ^ {mp/mehoVW^ix/xof^'-'+^^ » 1.(15) 

Although consideration of the post-shock fields is beyond the scope of the present 
paper, we can estimate the magnetic field spectrum at and after the shock jump as 
long as dissipation is not playing a role. The magnetic field of different correlation 
scales created in the foreshock is advected toward the shock, so a broad spectrum is 
accumulated: 

Ex oc A-S^ ~ A-° °9i, (16) 
where Eqs. ([9| and (10 1 were used and s = p — 1 ^ 1.2 was assumed. 



4. The afterglow foreshock 

The relation between the shock radius R and its Lorentz factor Fgh follows from a 
simple energy argument: the energy of an explosion is i? ~ (47r/3)P'^mpC^nisMrsh) 
therefore 

R ^ (10i« cm) eH^-^^'t;^'/' (17) 

or Fsh ~ ^52^'^isM^-^i8^^^i where E52 — E/IO^"^ erg and similarly for other quantities. 
The observed time of photons emitted at radius R is the time since the very first 
photons (i.e, emitted at P '--^ 0) arrive, tobs = R/c— R/vsh — R/c— P/[c(l — l/2rgj^)], 
that is ^obs ~ ^/(2rgijC). Using the equation for R, one gets 

r^h ^ 3.7 Eli\ii^%!^\ (18) 

R ^ (4.2 X 10^^ cm) Eli^-i^Xil. (19) 

Here we assumed a local GRB with z = 0; to include the redshift time dilation is 
trivial. 
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The co-moving density is no ^ nigMrgh (assuming the CR efficiency ^cR — 0-5 ~ 
1) and the minimum Lorentz factor of CR protons is 70 ^ Tgh- Hence, the length 
scale ~ Ao/(2^b) (Aq is the skin length) and the field Bq in the shock co-moving 
frame become 

xo ^{2x 10^ cm) n-^l^/{2(,B) - (10^ cm) n^^J.^ (20) 

Bo ^ (0.2 gauss) ^'J^nHl^T^i, ~ (1 gauss) E^'^R-.^/^ (21) 

where we assumed ^ 0.01. Assuming p — 2.2 and the interstellar fields, i?isM, to 
be of order a microgauss, we estimate X as 



X ^ xo 



0.2(^5 nisM)'/'i3rsL ~ 2 X 10^^ xo nfsM^fsM -6' (22) 



independent of T^h- On the other hand, 

i?/(2r,h) ~ (5 X 10^^ cm) Eli\ii^^T;^^/^ « X, (23) 

indicating that the ambient field is relatively unimportant, even for very steep energy 
spectra p ^ 3.5 rarely observed in prompt GRBs. The foreshock thickness is 

a^max - i?/(2r,h) ^ 5 X 10« Xo eI/,\-s]^^T;^^/\ (24) 
Therefore, the typical field in the foreshock is of sub-gauss strength: 

^(a^max) ~ (0.2 gauss) i^sV'^GM^rs (25) 
This field is relatively large-scale, as it's co- moving correlation scale is 

A(x„,ax) a;^ax/(2eB) - (5 X 10^^ cm) ^^s'/^'^iSM^if- (26) 
The power emitted by CR electrons from the foreshock amounts to 

P"*?^ ~ (2 X 10^9 erg s'^) E',^^LcR.45nTi^'Rli-° (27) 
in the observer's frame and is emitted at a peak (synchrotron) frequency 

- (10" Hz) £;|2°nrsM ^rs' (28) 
which corresponds to the IR band at about one day after the explosion, where R(t) is 



given in Eq. (|19|), so that cx tjay^^" ^day** 



5. Discussion 



Here we presented a model of a self-similar foreshock produced by protons scattered 
by a relativistic shock into the unmagnetized or weakly magnetized external medium. 
It is immediately applicable to the external shock producing a GRB afterglow. The 
model predicts that a large region in front of the shock, of thickness of order the 
shock radius, shall be filled with relatively strong and large-scale magnetic fields. 
The upstream magnetic field strength and correlation length depend on the distance 
from the shock and the power-law index of accelerated protons (cosmic rays) and are 
given by Eqs. ([9]), (10 1, (21) in the shock co-moving frame. The overall energetics 



of the field is dominated by large distances; hence the average foreshock B-field is of 
sub-milligauss strength and is increasing toward the shock while its typical coherence 
length is of order of few percent of the shock radius and is decreasing toward the 
shock. 

The result is interesting, especially in the light of observational constraints on 
the particle acceleration in GRB afterglows. It has been shown that a few-milligauss 
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magnetic fields are needed in front of the shock in order to efficiently Fermi-accelerate 



the electrons to the energies required to produce the observed X-ray emission (Li & 



Waxman 20061. Our model provides a possible and rather natural mechanism for 



generation of such fields in the far upstream medium. We also make a prediction 
that the shock-accelerated electrons will be radiating in the foreshock fields at a 
characteristic synchrotron frequency given by Eqs. (28 1. For reference, Um ~ 10 THz 
at about a day after the burst. It is possible that nonlinear effects omitted in this 
analysis (see discussion below) may limit the field strength to lower values compared 
to the present analysis, so the synchrotron peak can move to cm/mm-wave band. 
We can speculate that the emission from the foreshock can form an emission region 
separate from the afterglow shock and show up in the X-ray/optical band in the early 
afterglow phase and in the radio at the very late times. 

The presented results can have interesting implications for the radiative efficiency 
of external shocks. Unlike internal shocks, where electron cooling is extremely fast 
and a thin shock layer of thickness of a hundred ion inertial lengths may be enough 
to produce the observed prompt emission ( [Medvedev fc Spitkovsk;^|2008[ ) , the Weibel 
shock model (Medvedev & Loeb 19991 seems to face the efficiency problem when it 
is applied to an external shock. In such a shock, magnetic fields shall occupy a much 
larger region, perhaps the entire downstream region, in order for the shock to produce 
enough photons that will be observed as a delayed afterglow emission. This is not 
very likely (though not proved to be impossible, yet) provided that the small-scale 
fields generated at the shock can be subject to rapid dissipation. However, dissipation 
shall be of much lesser importance for the foreshock fields, which have a (much) larger 
coherence length, Eq. (|9|. In the steady state, the fields generated in the upstream 
are advected to the shock and their strength is maintained against dissipation by the 
anisotropy of the continuously "refreshed" CR distribution. Hence, one can expect 
that the magnetic field near the shock will have the spectrum given in Eq. ( 16 ). Once 
these fields pass though the shock into the downstream, they are enhanced by the 
shock compression and begin to decay. Commonly, dissipation is proportional to the 
inverse gradient scale squared oc 9^ cx A~^, so that the skin- length-scale fields may 
eventually disappear. However, the fields above a certain coherence length, Adiss, shall 
survive and fill up the post-shock medium. The mechanism of dissipation is not yet 
understood in detail, so it is premature to make any quantitative conclusions about 
A(jissj but it will certainly be much smaller than A(a;max) ^ R, see Eq. (26 1. Since 
Bx is a weak function of A one can speculate that relatively strong fields, perhaps of 
order tens or hundreds milligauss, to occupy the post-shock medium. 

We want to note that a number of simplifying assumptions has been made in 
our analysis. In particular, nonlinear feedback effects of the upstream magnetic field 
on the particle distribution, on the shock structure and on Fermi acceleration were 
omitted. The inclusion of these effects is hardly possible in any analytical model. 
We also assumed that a steady state exists for the shock-foreshock system at hand. 
Apparently, it is not at all clear whether the steady state is at all possible or the 
system exhibits an intermittent behavior. One can envision a scenario in which the 
CRs overproduce upstream magnetic fields leading to enhanced particle scattering and 
the overall preheating of the ambient medium, which, in turn, can cause the shock to 
weaken, disappear and then re-appear in a different place further upstream. Presently 
available 2D PIC simulations of an electron-position shock do show the upstream field 
amplification and no steady state has been achieved: both upstream and downstream 
fields continue to grow for the duration of the simulations ( Keshet et al. 2008 1 . We 
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argue that extensive PIC or/and hybrid simulations of a shock are imperative for 
further study. 

Finally, we mention that our model compliments other models of the magnetic 
field generation. It is reasonable to expect that the field can be amplified by vortical 
motions produced by the Richtmeyer-Meshkov instability, if the ambient medium is 



clumpy or if the shock velocity is not perfectly uniform ( Goodman & MacFadyen 2007 



magnetic fields are stron; 


I enough, the fields can be generated via non-resonant (Bell 


& Lucek||2001 


Bell 


2004 


Pelletier et al.||2008l or resonant (Diamond & Malkov||2006 


2007 Zweibel 


2003 


Vladimirov et al.||2006 Ellison et al.||2007l mechanisms, or both 


(Reville et al. 


2008 
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Figure 1. A schematic representation of the foreshock magnetic fields: the 
coherence length is increasing with the upstream distance. Below are schematic 
graphs showing variation of the spectrum of the streaming part of cosmic rays 
and the corresponding self-generated fields (highlighted). 



